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Abstract 



We have generalized the approach of Brodsky et al. for the intrinsic charm quark distribution in the nucleons to the 
light-quark sector involving intrinsic u, d, s and s sea quarks. We compare the calculations with the existing d — u, s + s, 
CNI and u + d — s — s data. The good agreement between the theory and the data allows the extraction of the probabilities 
for the \uuduu), \uuddd), and \uudss) five-quark Fock states in the proton. We also calculate the ^-dependence of the 
intrinsic charm after taking into consideration the QCD evolution of the intrinsic quark distribution. 



o 



CO 



Key words: nucleon structure, five-quark states, intrinsic sea quark 
PACS: 12.38.Lg, 14.20.Dh, 14.65.Bt, 14.65.Dw, 13.60.Hb 



The origin of sea quarks of the nucleons remains a 
1 subject of intense interest in hadron physics. Brodsky, 



'Q^Hoyer, Peterson, and Sakai (BHPS) [1| suggested some 
I time ago that there are two distinct components of the 
nucleon sea. The first is called the "extrinsic" sea origi- 
pH nating from the splitting of gluons into QQ pairs. This 
i— i extrinsic sea can be well described by quantum chromo- 
dynamics (QCD). Another component of the nucleon sea 
^ is the "intrinsic" sea which has a nonperturbative origin. 

I In particular, the \uudQQ) five-quark Fock states can lead 
00 "to the "valence-like" intrinsic sea for the Q and Q in the 
CO proton. This intrinsic component is expected to carry a 
relatively large momentum fraction x, in contrast to the 
extrinsic one peaking at the small- a: region. Brodsky et 
'al. [H proposed that the \uudcc) five-quark state can lead 
| | !to enhanced production of charmed hadrons at the forward 
. . rapidity region. The CTEQ collaboration [2j has exam- 
^ ined all relevant hard-scattering data and concluded that 
the data are consistent with a wide range of the intrinsic 
.charm magnitude, ranging from null to 2-3 times larger 
Cythan the estimate by the BHPS model. This suggests that 
more precise experimental measurements are needed for 
determining the magnitude of the intrinsic charm compo- 
nent. 

In a recent work Q, we generalized the BHPS model of 
the five-quark Fock states to the light-quark sector. This 
work was motivated by the expectation that the proba- 
bility for the \uudQQ) Fock state is approximately pro- 
portional to l/rrig, where itiq is the mass of the quark 
Q Although this 1/toq dependence is applicable only 
when the quark mass is heavy the light five-quark 
states \uuduu), \uuddd) and \uudss) are likely to have sig- 
nificantly larger probabilities than the \uudcc) state, and 
could be more readily observed experimentally. 

By solving the Bjorken-x distribution of the Q sea 



quark for the \uudQQ) five-quark state in the BHPS model 
numerically, it was found [3[ that the existing d(x) — u(x) 
and u{x) + d(x) — s(x) — s(x) data can be well described 
by the calculation, provided that the QCD evolution [j| 
of these distributions is taken into account. Moreover, the 
probabilities for the \uuduu) and the \uuddd) five-quark 
states could also be extracted from these data. However, 
the extracted values of these two probabilities depend on 
the assumption adopted for the probability of the \uudss) 
state [H. 

In this paper, the previous work is extended further to 
determine the probability of the \uudss) five-quark state 
using the recent s{x) + s(x) data from the HERMES col- 
laboration Q. We found that the s(x) + s(x) data in 
the x > 0.1 region are quite well described by the BHPS 
model, allowing the extraction of the probability of the 
\uudss) state. Using this probability for the \uudss) five- 
quark component, more precise values for the \uuduu) and 
the \uuddd) states could then be obtained from the com- 
parison of the BHPS calculations with the d(x) —u(x) and 
u(x) + d(x) — s(x) — s(x) data. We have also examined 
the effect of the QCD evolution on the x distribution of 
the intrinsic charm. In particular, we note that the region 
most sensitive to intrinsic charm is shifted to lower x as a 
result of QCD evolution. This has implication on future 
searches for intrinsic charm. 

For a \uudQQ) five-quark Fock state of the proton, the 
probability for quark i to carry a momentum fraction Xi 
is given in the BHPS model [1] as 



P(x u ...,x 5 ) = N 5 5(l 
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where the delta function ensures that the proton momen- 
tum is shared among the individual constituents. A5 is 
the normalization factor for the five-quark Fock state, and 
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rrii is the mass of quark i. Eq. [T]was solved analytically in 
Ref. [l[ for the limiting case of 7714,5 >> m pi m i,2,3, where 
m p is the proton mass. For the more general case, Eq. [T] 
can be solved numerically as discussed in Ref. Q. In par- 
ticular, the x distribution of Q in the \uudQQ) state, called 
P®Q(xq), can be calculated numerically. The moment of 

PQQ(xq) is defined as V®® , namely, 



V, 



QQ 



P^{x Q )dx Q . 



(2) 



V®® represents the probability of the \uudQQ) five-quark 
Fock state in the proton. In the limit of 7714 5 >> m p , 7711,2,3, 
one can obtain = V 5 /(3600m| 5 ). For the more 

general case, the relation between V®® and N$ can be 
calculated numerically 

To compare the experimental data with the prediction 
based on the intrinsic five-quark Fock state, it is necessary 
to separate the contributions of the intrinsic sea quark and 
the extrinsic one. The d(x) — u(x) is an example of quan- 
tities which are free from the contributions of the extrinsic 
sea quarks, since the perturbative g — > QQ processes will 
generate uu and dd pairs with equal probabilities and have 
no contribution to this quantity. The d(x)—u(x) data from 
the Fermilab E866 Drell-Yan experiment at the Q 2 scale of 
54 GeV 2 [7| are shown in Fig. [T] Also shown in Fig. [T] are 
the data obtained at a lower scale of Q 2 = 2.5 GeV 2 by the 
HERMES collaboration in a semi-inclusive deep-inelastic 
scattering (SIDIS) experiment [8j. 

The BHPS model has a specific prediction on the shapes 
of the x distributions for d and u, since these anti-quarks 
originate from the | uuddd) and | uuduu) configurations and 
can be readily calculated. In the BHPS model, the u and d 
are predicted to have the same x-dependence if m u — rrid- 
However, the probabilities of the \uuddd) and \uuduu) con- 
figurations, 7-5 and P5", are not known from the BHPS 
model, and remain to be determined by the experiments. 
Non-perturbative effects such as Pauli-blocking IS] could 
lead to different probabilities for the \uuddd) and \uuduu) 
configurations. Nevertheless the shape of the d(x) — u(x) 
distribution shall be identical to those of d{x) and u(x) in 
the BHPS model. Moreover, the normalization of d(x) — 
u(x) is known from the measurement of Fermilab E866 
Drell-Yan experiment Q as 



{d{x)-u{x))dx = V, 



dd 
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0.118±0.012.(3) 



Equation [3] allows us to compare the calculations from the 
BHPS model with the d(x) — u(x) data. 

The d{x) — u{x) distribution from the BHPS model is 
first calculated using Eq. [T]with m u = 777^ = 0.3 GeV/c 2 , 
and m p = 0.938 GeV/c 2 , and Eq. [3]for the normalization. 
Since the E866 and the HERMES data were obtained at 
Q 2 of 54 GeV 2 and 2.5 GeV 2 , respectively, it is important 
to evolve the d(x) —u{x) distribution from the initial scale 
/j, expected to be around the confinement scale, to the Q 2 
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Figure 1: Comparison of the d(x)—u(x) data from Fermilab E866 and 
HERMES with the calculations based on the BHPS model. Eq. Q] 
and Eq. [3] were used to calculate the d(x) — u(x) distribution at 
the initial scale. The distribution was then evolved to the Q 2 of 
the experiments and shown as various curves. Two different initial 
scales, fi = 0.5 and 0.3 GeV, were used for the E866 calculations in 
order to illustrate the dependence on the choice of the initial scale. 



corresponding to the data. As d[x) — u(x) is a flavor non- 
singlet parton distribution, its evolution from /1 to Q only 
depends on the values of d(x) — u(x) at the initial scale /j, 
and can be readily calculated using the non-singlet evolu- 
tion equation [5j. For the initial scale, we adopt the value 
of n = 0.5 GeV, which was chosen by Gliick, Reya, and 
Vogt [10( in the so-called "dynamical approach" using only 
valence-like distributions at the initial fj, 2 scale and relying 
on evolution to generate the quark and gluon distributions 
at higher Q 2 . 

The solid and dashed curves in Fig. [1] correspond to 
d{x) — u(x) calculated from the BHPS model evolved to 
Q 2 = 54 GeV 2 using fi = 0.5 and 0.3 GeV, respectively. 
The x-dependence of the E866 d(x)—u(x) data is quite well 
described by the five-quark Fock states in the BHPS model 
provided that the Q 2 -evolution is taken into consideration. 
It is interesting to note that an excellent fit to the data 
can be obtained if /1 = 0.3 GeV is chosen (dashed curve in 
Fig.[TJ rather than the more conventional value of /1 = 0.5 
GeV. Also shown in Fig. [T] are the calculations with the 
BHPS model evolved to Q 2 = 2.5 GeV 2 using /1 = 0.5 GeV 
and /i = 0.3 GeV. The calculations are in agreement with 
the HERMES data within the experimental uncertainties. 

We now consider the extraction of the \uudss) five- 
quark component from existing data. The HERMES col- 
laboration reported the determination of x(s(x) + s(x)) 
over the range of 0.02 < x < 0.5 at Q 2 = 2.5 GeV 2 from 
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Figure 2: Comparison of the HERMES x(s(x) + s(x)) data with the 
calculations based on the BHPS model. The solid and dashed curves 
are obtained by evolving the BHPS result to Q 2 = 2.5 GeV 2 using 
/i = 0.5 GeV and fi = 0.3 GeV, respectively. The normalizations of 
the calculations are adjusted to fit the data at x > 0.1 with statistical 
errors only, denoted by solid circles. 

their measurement of charged kaon production in SIDIS re- 
action [6]. The HERMES data, shown in Fig. EJ exhibits 
an intriguing feature. A rapid fall-off of the strange sea 
is observed as x increases up to x ~ 0.1, above which the 
data become relatively independent of x. The data suggest 
the presence of two different components of the strange 
sea, one of which dominates at small x (x < 0.1) and the 
other at larger x (x > 0.1). This feature is consistent 
with the expectation that the strange-quark sea consists 
of both the intrinsic and the extrinsic components hav- 
ing dominant contributions at large and small x regions, 
respectively. In Fig. 2 we compare the data with calcula- 
tions using the BHPS model with m s = 0.5 GeV/c 2 . The 
solid and dashed curves are results of the BHPS model 
calculations evolved to Q 2 = 2.5 GeV 2 using /i — 0.5 GeV 
and fj, = 0.3 GeV, respectively. The normalizations are 
obtained by fitting only data with x > 0.1 (solid circles in 
Fig. E]) , following the assumption that the extrinsic sea has 
negligible contribution relative to the intrinsic sea in the 
valence region. Figure 2 shows that the fits to the data are 
quite adequate, allowing the extraction of the probability 
of the \uudss) state as 

Vf = 0.024 {fn = 0.5 GeV); 

Pf = 0.029 (/i = 0.3 GeV). (4) 

We consider next the quantity u(x) + d(x) — s(x) — 
s(x). Combining the HERMES data on x(s(x) + s(x)) with 



Figure 3: Comparison of the x(d(x)+u(x) — s(x) — s(x)) data with the 
calculations based on the BHPS model. The values of x(s(x) + s(x)) 
are from the HERMES experiment [(J, and those of x(d(x) + ii(x)) 
are obtained from the PDF set CTEQ6.6 [lj. The solid and dashed 
curves are obtained by evolving the BHPS result to Q 2 = 2.5 GeV 2 
using fA = 0.5 GeV and /i = 0.3 GeV, respectively. The normalization 
of the calculations are adjusted to fit the data. 

the x(d(x) + u(x)) distributions determined by the CTEQ 
group (CTEQ6.6) pJ| . the quantity x(u(x) + d(x) — s(x) - 
s(x)) can be obtained and is shown in Fig. [21 This ap- 
proach for determining x(u(x) + d(x) — s(x) — s(x)) is iden- 
tical to that used by Chen, Cao, and Signal in their recent 
study of strange quark sea in the meson-cloud model [l2| • 
An important property of u + d — s — s is that the 
contribution from the extrinsic sea vanishes, just like the 
case for d~u. Therefore, this quantity is only sensitive to 
the intrinsic sea and can be compared with the calculation 
of the intrinsic sea in the BHPS model. We have 

u(x) + d(x) — s(x) — s(x) — 
P ua (x a ) + P dct (x s ) - ■>!'"{.,■ J. (5) 

We can now compare the x(u{x) + d(x) — s(x) — s(x)) 
data with the calculation using the BHPS model. Since 
u + d— s — sis a flavor non-singlet quantity, we can readily 
evolve the BHPS prediction to Q 2 = 2.5 GeV 2 using fx = 
0.5 GeV and the result is shown as the solid curve in Fig.[3J 
It is interesting to note that a better fit to the data can 
again be obtained with fj, = 0.3 GeV, shown as the dashed 
curve in Fig. [3] 

From the comparison between the data and the BHPS 
calculations shown in Figs.[T][3J we can determine the prob- 
abilities for the \uuduu), \uuddd), and \uudss) configura- 
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tions as follows: 



VT = 0.162; 



Vf = 0.240; 



Vf = 0.280; 



7> 5 SS = 0.024 
(H = 0.5 GeV) 



VI s = 0.029 
(fj. = 0.3 GeV) 



(6) 



(7) 



l O 



depending on the value of the initial scale fi. It is re- 
markable that the d(x) — u(x), the s(x) + s(x), and the 
d[x) + u(x) — s(x) — s(x) data not only allow us to check 
the predicted x-dependence of the five-quark Fock states, 
but also provide a determination of the probabilities for 
these states. 

Equations 6 shows that the combined probability for 
proton to be in the \uudQQ) states is around 40%. It is 
worth noting that an earlier analysis of the d—u data in the 
meson cloud model concluded that proton has ~60%_prob- 



ability to be in the three-quark bare-nucleon state 
qualitative agreement with the finding of this study. A sig- 
nificant feature of the present work is the extraction of the 
\uudss) component, which would be related to the kaon- 
hyperon states in the meson cloud model. It is also worth 
mentioning that in the BHPS model the \uudQQ) states 
have the same contribution to the proton's magnetic mo- 
ment as the \uud) three-quark state, since Q and Q in the 
\uudQQ) states have no net magnetic moment. Therefore, 
the good description of the nucleon's magnetic moment 
by the constituent quark model is preserved even with the 
inclusion of a sizable five-quark components in the BHPS 
model. 

We note that the probability for the \uudss) state is 
smaller than those of the \uuduu) and the \uuddd) states. 
This is consistent with the expectation that the probability 
for the \uudQQ) five-quark state is roughly proportional 
to 1 /rriQ [l|, |4[ . One can then estimate that the probability 
for the intrinsic charm from the \uudcc) Fock state, V^ c to 
be roughly 0.01. This is also consistent with an estimate 



based on the bag model [14J, as well as with an analysis 
of the EMC charm-production data [lij]. Figure 4 shows 
the x distribution of intrinsic c calculated with the BHPS 
model using 1.5 GeV/c 2 for the mass of the charm quark. 
Also shown in Fig. 4 is the calculation which evolve the 
BHPS calculation from the initial scale, n = 0.5 GeV, to 
Q 2 = 75 GeV 2 , the largest Q 2 scale reached by EMC [L$. 
It is interesting to note that the intrinsic charm contents 
at the large x (x > 0.3) region are drastically reduced 
when Q 2 evolution is taken into account. Figure 4 suggests 
that the most promising region to search for evidence of 
intrinsic charm could be at the somewhat lower x region 
(0.1 < x < 0.4), rather than the largest x region explored 
by previous experiments. It is worth noting that we adopt 
the simple assumption that the initial scale is the same for 
all five-quark states. It is conceivable that the initial scale 
for intrinsic charm is significantly higher due to the larger 
mass of the charmed quark. The dashed curve shows the x 
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Figure 4: Calculations of the c(x) distributions based on the BHPS 
model. The solid curve corresponds to the calculation using Eq. [l] 
and the dashed and dotted curves are obtained by evolving the BHPS 
result to Q 2 = 75 GeV 2 using fi = 3.0 GeV, and fi = 0.5 GeV, 
respectively. The normalization is set at VS C = 0.01. 



distribution of intrinsic c at Q 2 = 75 GeV 2 when the initial 
scale is set at fi = 3 GeV, corresponding to the threshold 
of producing a pair of charmed quarks. As expected, the 
shape of the intrinsic c x distribution becomes similar to 
that of the BHPS model. 

In conclusion, we have generalized the existing BHPS 
model to the light-quark sector and compared the calcu- 
lation with the d — u, s + s, and u + d — s — s data. The 
qualitative agreement between the data and the calcula- 
tions provides strong support for the existence of the in- 
trinsic u, d and s quark sea and the adequacy of the BHPS 
model. This analysis also led to the determination of the 
probabilities for the five-quark Fock states for the proton 
involving light quarks only. This result could guide future 
experimental searches for the intrinsic c quark sea or even 
the intrinsic b quark sea (l7| . which could be relevant for 
the production of Higgs boson at LHC energies [l8[ . 
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